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LONG-TERM  GOALS 

The  over-all  goal  of  this  research  is  the  development  of  an  accurate  and  reliable  propagation  model 
applicable  to  environments  which  exhibit  strong  range  dependence  in  all  three  spatial  dimensions. 

OBJECTIVES 

The  objective  of  this  work  is  to  gain  an  understanding  the  physics  of  propagation  in  continental  shelf 
areas,  specifically  horizontal  refraction  and  mode  coupling  induced  by  three-dimensional  (3D) 
inhomogeneities  in  the  waveguide.  A  coupled-mode  approach  has  been  applied  for  this  purpose.  The 
coupled-mode  approach  is  attractive  for  solving  problems  involving  3D  propagation  for  several  reasons. 
First,  this  technique  provides  intuitive  results  for  understanding  the  features  responsible  for  observed 
propagation  effects  in  range-dependent  environments.  For  example,  upslope  propagation  is 
characterized  by  acoustic  energy  radiated  into  the  bottom  at  discrete  depths  associated  with  mode 
cut-off.  The  modal  decomposition  of  the  acoustic  field  has  also  been  used  to  describe  horizontal 
refraction  in  a  wedge-shaped  ocean,  for  which  the  single-mode  interference  pattern  associated  with  rays 
launched  up  and  across  the  shelf  has  been  well  documented  [Weinberg  and  Burridge  (1974)]. 
Furthermore,  coupled-mode  solutions  are  highly  accurate  and  have  been  used  for  benchmarking 
solutions  to  range-dependent  problems  [Jensen  and  Ferla  (1990)].  In  order  to  appreciate  the  limitations 
of  existing  3D  models,  it  is  necessary  to  have  methods  which  can  provide  reference  solutions  for 
comparison. 

APPROACH 

A  3D  acoustic  propagation  model  based  on  the  stepwise  coupled-mode  approach  [Evans  (1983)] 
implemented  as  a  single-scatter  solution  has  been  developed.  This  technique  is  based  on  a  hybrid 
modeling  approach  for  which  normal  modes  are  applied  in  the  vertical  dimension  and  a  PE  solution  is 
applied  in  the  horizontal  dimension. 

The  inhomogeneous  Helmholtz  equation  for  pressure  P(r.  O.z.)  at  range  r,  azimuth  0,  and  depth  z  from  a 
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point  continuous  wave  source  of  amplitude  S(co)  located  at  range  r  =  0  and  depth  z  —  zo,  is  given  by 


p(r,0,z)V. 


lP(r,9,z) 


VP{r,e,z) 


+  k2(r ,  6,z)P(r,  d,z )  =  -47rS(©)^— ^-<5(z-zo). 


(1) 


where  k(r.  O.z)  —  Co/c(r ,  O.z)  is  the  acoustic  wavenumber,  (0  =  2k  f,  f  is  the  acoustic  frequency, 
c(r,  0,z)  is  sound  speed,  and  p(r,  0,z)  is  density. 

The  solution  for  pressure  is  found  by  a  separation  of  variables 


M 

P(r,  0,z)  =  0)<l>m(z;r,  0), 

777=1 


(2) 


where  Am(r,  0)  are  the  modal  amplitudes,  and  (j)m  (z:,  r.  0)  are  the  modal  eigenfunctions. 
The  modal  eigenfunctions  (j)m(z.',r,  0)  satisfy 


p(r,0,z)  — 
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p(r,0,z)  dz 


+  [£2 0,  9 ,  z)  -  k;n  (r,  0 )]  <j)m  (z;  r,  0 )  =  0, 


(3) 


where  km(r,  0)  is  the  horizontal  wavenumber  of  the  in11'  mode.  Boundary  conditions  are  defined  by  the 
plane  wave  reflection  coefficient  below  the  upper  halfspace  and  above  the  lower  halfspace.  In  this  work, 
the  Pekeris  branch  cut  is  chosen  such  that  the  total  pressure  is  calculated  from  a  Pekeris  branch  line 
integral,  plus  a  finite  sum  of  trapped  modes,  plus  an  infinite  sum  of  leaky  modes.  A  small  gradient  is 
introduced  in  the  lower  halfspace  which  effectively  removes  the  branch  point  and  associated  branch  cut 
from  the  problem  [Westwood  and  Koch  (1999)].  As  a  result,  the  leaky  modes  eventually  decay  as  a 
function  of  depth  in  the  lower  half  space.  The  eigenfunctions  are  normalized  so  that 
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In  the  absence  of  mode  coupling,  the  adiabatic  solution  for  the  mode  amplitudes  Am(x,y)  satisfies 
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This  equation  must  be  solved  for  each  mode  with  the  horizontal  refraction  determined  by  the  modal 
phase  speed  cpi,m (x, y)  =  c o/Re{k,n(x,y )}  and  modal  attenuation  am (x, y)  =  Im {km(x,y)}.  Such  a  two 
dimensional  (2D)  Helmholtz  equation  can  be  solved  by  standard  techniques.  In  this  work,  the  solution 
to  the  horizontal  refraction  equation  (Eq.  (5))  is  obtained  from  a  PE  model  [Collins  (1994)]  which  has 
been  modified  for  cylindrical  coordinates  [Napolitano  (1985)], 


^7  =  ik0  \/l+k02(^2j^2+k2-ki)Am,  (6) 

where  ko  —  ®  is  the  reference  wavenumber,  and  cq  is  the  reference  sound  speed. 

Mode-coupling  is  incorporated  into  the  PE  solution  by  a  stepwise  coupled-mode  technique  [Evans 
(1983)].  This  approach,  originally  derived  for  the  2D  Helmholtz  equation,  is  applied  by  discretizing  a 
range-dependent  environment  into  a  series  of  range-independent  segments  for  which  continuity  of 
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pressure  and  particle  velocity  must  be  satisfied  at  the  vertical  boundaries  between  segments. 
Application  of  the  boundary  conditions  results  in  a  set  of  linear  equations  from  which  the  coupling 
matrix  Rm,n(r ,  9)  is  obtained, 


RmJr-  9)  =  C*Ln(r,  9) +C^Kn(r,  0) 


•LR 
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where  the  suffixes  L  and  R  denote  the  properties  to  the  left  and  right  of  a  vertical  interface.  This  form  of 
the  mode-coupling  matrix  describes  the  single-scatter  approximation,  such  that  each  pair  of  segments  is 
treated  as  an  individual  problem,  thus  neglecting  higher-order  terms  resulting  from  multiple  scattering 
at  other  interfaces.  It  can  be  shown  that  in  the  absence  of  mode  coupling,  the  mode-coupling  matrix 
reduces  to  the  identity  matrix,  and  the  adiabatic  solution  is  obtained. 


In  the  3D  coupled-mode  model,  the  range-independent  segments  are  defined  in  the  radial  direction, 
forming  angular  sectors  in  the  range -bearing  plane.  By  incorporating  the  mode-coupling  matrix  into  Eq. 
(6),  the  solution  for  the  coupled-modal  amplitudes  is  obtained, 


dAm 

dr 


N  \  d2 

Rm,n{p  9 )An  -\-  ik0  1  T ~^2  d~k~  k q  Am. 


(9) 


According  to  this  solution  technique,  at  each  range  step,  mode-coupling  is  calculated  by  the  first  term  in 
Eq.  (9)  and  horizontal  refraction  is  calculated  by  the  second  term.  Although  mode-coupling  occurs  only 
in  the  radial  direction,  the  coupled  energy  is  refracted  out  of  the  range-depth  plane  by  the  second  term. 


WORK  COMPLETED 

The  main  accomplishments  of  2012  include:  (1)  a  comparison  of  the  3D  coupled-mode  solution  to  a 
reference  solution  obtained  with  semi-analytic  technique  for  propagation  around  a  shallow-water 
seamount,  (2)  the  addition  of  a  3D  rough  sea  surface  in  the  model  with  comparison  to  solution 
calculated  with  a  3D  PE  model,  and  (3)  an  implementation  of  parallel  processing  to  allow  for  efficient 
calculations  at  higher  frequencies  and  longer  ranges. 


RESULTS 

Propagation  around  a  conical  seamount 

The  3D  coupled-mode  model  described  above  is  applied  to  calculate  propagation  around  a  conical 
seamount  in  shallow  water.  This  environment  was  chosen  because  the  shallow-water  seamount  induces 
both  strong  horizontal  refraction  and  mode-coupling  effects.  Additionally,  a  semi-analytic  solution  is 
available  for  comparison  [Luo  and  Schmidt  (2009)].  The  waveguide  is  250  m  deep,  bounded  above  by  a 
flat  sea  surface  and  below  by  a  penetrable  bottom.  The  water  column  sound  speed,  density,  and 
attenuation  are  1500  m/s,  1.0  g/cm3,  and  0.0  dB/A.  The  seamount  rises  to  a  height  of  100  m  above  the 
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(a)  3D  Coupled  Mode 


(b)  3D  Adiabatic  Mode 


(c)  Nx2D  Coupled  Mode 
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Figure  1:  ( Color  online)  Transmission  loss  at  a  depth  of  100  m  calculated  by  3D  (a)  coupled-  and  (b) 
adiabatic-mode  techniques  and  by  Nx2D  (c)  coupled-  and  (d)  adiabatic -mode  techniques. 


seafloor  and  has  a  radius  of  350  m  at  its  base.  The  seabed  sound  speed,  density,  and  attenuation  are 
1800  m/s,  2.0  g/cm3,  and  0.1  dB/A.  The  halfspace  includes  an  attenuation  gradient  such  that 
Aa  =  ayA  —  cchs  =  0.01  dB/A,  where  <Xhs  is  the  attenuation  of  the  halfspace,  and  a:y  is  the  attenuation 
three  wavelengths  into  the  halfspace.  The  acoustic  source  is  located  800  m  from  the  axis  of  the 
seamount  at  a  depth  100  m.  The  solution  was  calculated  for  an  excitation  frequency  of  40  Hz.  At  this 
frequency,  there  are  seven  trapped  modes  at  the  location  of  the  source.  Adding  five  leaky  modes  was 
sufficient  to  obtain  a  convergent  solution. 

Transmission  loss  at  a  depth  of  100  m  calculated  by  four  different  approximations  is  shown  in  Figure  1. 
The  effects  of  horizontal  refraction  are  assessed  by  comparing  results  from  3D  and  Nx2D  models, 
which  neglect  out-of-plane  effects  by  applying  a  2D  model  to  predict  TL  for  vertical  slices  through  the 
environment  [Perkins  and  Baer  (1982)].  The  effects  of  mode-coupling  are  assessed  by  comparing 
adiabatic-  and  coupled-mode  solutions. 

The  importance  of  mode-coupling  is  observed  from  the  differences  in  the  coupled-  and  adiabatic-mode 
solutions  shown  in  Figure  1(a)  and  (b).  In  these  figures,  there  are  significant  differences  in  the  acoustic 
field  in  the  region  on  the  far  side  of  the  seamount.  The  effects  of  horizontal  refraction  can  be  understood 
by  comparing  the  3D  to  the  Nx2D  solutions.  For  the  coupled-mode  solutions  (Figure  1(a)  and  (c)), 
differences  in  the  modal  interference  patterns  are  observed  behind  the  seamount.  In  addition,  the  3D 
solution  is  characterized  by  a  wider  diffractive  pattern  behind  the  seamount,  as  sound  is  refracted  away 
from  the  seamount  and  into  deeper  water.  The  effects  of  horizontal  refraction  can  be  viewed  even  more 
clearly  in  the  adiabatic-mode  solutions  (Figure  1(b)  and  (d)).  For  this  case,  the  3D  solution  shows 
higher  loss  directly  behind  the  seamount,  evidence  of  sound  being  refracted  away  from  the  seamount. 
Finally,  both  3D  solutions  show  fringes  on  the  outer  edges  of  the  affected  area.  These  result  from 
single-mode  interference  patterns  between  refracted  and  non-refracted  acoustic  paths. 
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To  evaluate  the  accuracy  of  the  3D  coupled-mode  solution,  it  is  compared  with  a  published  result 
calculated  using  a  semi-analytic  technique  [Luo  and  Schmidt  (2009)].  According  to  this  method,  the 
conical  seamount  is  divided  into  a  number  of  rings,  in  each  of  which  a  series  expansion  of  the  acoustic 
pressure  in  terms  of  normal  modes  and  an  azimuthal  Fourier  series  is  applied.  The  3D  coupled-mode 
solution  (Figure  2(a))  shows  excellent  quantitative  agreement  with  the  published  result  (see  Figure  6(c) 
of  [Luo  and  Schmidt  (2009)]).  The  features  described  for  the  3D  coupled-mode  solution  can  also  be 
observed  in  the  semi-analytic  solution  by  the  wider  shadow  behind  the  seamount  and  by  the  interference 
fringes  flanking  the  affected  region. 

Propagation  under  a  rough  sea  surface 

To  demonstrate  the  capability  of  the  3D  coupled-mode  model  to  calculate  TL  under  a  rough  sea  surface, 
an  environment  with  a  2D  sinusoidal  sea  surface  is  considered.  Scattering  from  2D  sinusoidal  surfaces 
has  been  studied  extensively  for  the  2D  propagation  problem,  for  which  the  surface  is  composed  of 
curved,  axisymmetric  wavefronts,  and  more  recently  for  the  3D  propagation  problem,  for  which  the 
surface  is  composed  of  planar  wavefronts.  The  solution  from  the  3D  coupled-mode  model  is  compared 
to  results  calculated  by  a  3D  wide  angle  PE  (WAPE)  model  [Smith  (2012)]. 

For  the  problem  considered  in  this  work,  the  height  of  the  sea  surface  is  described  by 

s(x)  =Asin(-^^  +  y)  (10) 

A'SU.rf  z- 

where  A  =  10  m  is  in  the  amplitude  of  the  surface  fluctuations,  and  Xsurf  =  200  m  is  the  wavelength  of 
the  fluctuations.  The  surface  fluctuations  form  corrugations  in  the  horizontal  plane  which  are 
independent  of  the  y-direction.  The  mean  water  depth  is  100  m.  The  water  column  sound  speed, 
density,  and  attenuation  are  1500  m/s,  1.0  g/cm3,  and  0.0  dB/A.  The  seabed  sound  speed,  density,  and 
attenuation  are  1700  m/s,  1.5  g/cm3,  and  0.15  dB/A.  The  halfspace  includes  an  attenuation  gradient 
such  that  A  a  =  a^x  —  <Xhs  =  0.01.  Transmission  loss  was  calculated  for  a  1  kHz  source  located  at  30  m 
water  depth.  The  solution  was  calculated  from  a  sum  of  99  modes.  Including  additional  higher-order 
modes  in  the  solution  (up  to  199  modes)  was  found  to  have  no  significant  effect  on  the  solution  beyond 
a  range  of  0.1  km. 

As  shown  in  Figure  2,  there  are  significant  differences  in  the  3D  coupled-  and  adiabatic-mode  solutions, 
indicating  the  importance  of  mode-coupling  effects  in  the  solution.  As  observed  from  Figure  2(a)  the 
coupled-mode  solution  shows  vertical  refraction  of  sound  under  the  first  wave  crest,  whereas  the 
adiabatic-mode  solution  (Figure  2(c))  fails  to  capture  scattering  in  the  vertical  plane.  Further  differences 
are  apparent  in  the  range-bearing  plots  (Figure  2(b,d))  where  the  interference  patterns  differ,  especially 
around  x  =  ±0.2  km  where  the  vertical  scattering  of  sound  under  the  first  wave  trough  occurs  in  the 
coupled-mode  solution. 

The  3D  coupled-mode  solution  (Figure  2(a,b))  shows  many  of  the  same  features  as  the  solution  from  the 
3D  WAPE  approximation  (see  Figure  2  of  [Smith  (2012)]).  For  example,  both  the  coupled-mode  and 
PE  solutions  are  able  to  capture  the  vertical  scattering  caused  by  interaction  with  the  rough  sea  surface. 
Additionally,  in  the  range-bearing  slice  at  30  m  depth,  both  solutions  show  interference  patterns  which 
form  oblong  rings  near  the  source  and  transition  to  parabolic-shaped  features  for  distances  greater  than 
y  =  0. 1  km.  The  aforementioned  line  of  interference  patterns  centered  around  around  x  =  ±0.2  km  is 
also  visible  in  both  solutions.  However,  there  are  some  significant  differences  in  the  solutions  obtained 
by  the  two  techniques.  The  result  calculated  by  the  3D  WAPE  model  predicts  significantly  lower  TL 
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Figure  2:  ( Color  online)  Transmission  loss  under  the  sinusoidal  sea  surface  from  (a,b)  coupled-  and 
(c,d)  adiabatic-mode  solutions:  (a,c)  range-depth  slice  ofTL  along  single  radial  ( perpendicular  to 
surface  wave  displacement);  (b,d)  range-bearing  slice  at  30  m  depth. 


than  the  Nx2D  solution  (see  Figure  1  of  [Smith  (2012)])  in  the  direction  perpendicular  to  the  wave 
crests.  On  the  other  hand,  the  3D  coupled-mode  solution  contains  small  differences  (less  than  1  dB  at 
most  depths  and  ranges)  from  the  Nx2D  solution  along  the  x-axis.  These  differences  are  evident  in  the 
range-bearing  dependence  as  well,  especially  for  ranges  greater  than  0.2  km. 

High  Performance  Computing  Capabilities 

Coupled-mode  codes  are  intrinsically  computationally  intensive.  Efficient  computation  of  the  solution 
was  achieved  by  parallel  modal  computations,  which  are  performed  simultaneously  over  azimuth  at 
each  range  step.  Message  Passing  Interface  (MPI)  was  implemented  in  the  code  for  this  purpose. 
Computation  time  is  proportional  to  the  number  of  processors  utilized  by  the  computation.  The  result  is 
a  significant  reduction  in  computation  time,  such  that  parallelized  computations  can  be  completed 
hundreds  of  times  faster  than  serial  computations.  This  has  made  it  possible  to  compute  solutions 
problems  involving  higher  frequencies  and  longer  ranges.  The  parallelized  code  is  highly  portable  and 
has  been  run  on  multiple  platforms  including  the  Sun  Constellation  Linux  Cluster  Ranger  [The 
University  of  Texas  at  Austin  (2012)]  and  the  U.S.  Air  Force  Research  Laboratory  (AFRL)  DoD 
Supercomputing  Resource  Center  (DSRC)  [U.S.  Air  Force  Research  Laboratory  (2012)]. 

IMPACT/APPLICATIONS 

The  impact  of  this  work  will  be  an  increased  understanding  of  acoustic  propagation  through 
complicated  coastal  environments  for  which  the  bathymetry,  seabed  properties,  and  oceanography  can 
vary  in  three  dimensions. 
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TRANSITIONS 


The  primary  transition  for  this  project  is  an  accurate  and  reliable  model  for  acoustic  propagation  in 
environments  with  strong  three-dimensional  range  dependence.  Because  coupled-mode  approaches  are 
computationally  intensive,  they  have  historically  been  used  to  benchmark  faster  techniques  which 
approximate  the  solution  to  the  wave  equation. 

RELATED  PROJECTS 

Geoacoutic  inversion  in  three-dimensional  environments 

The  goal  of  this  project  is  to  estimate  water  column  sound  speed  in  a  3D  volume  using  modal  travel 
time  measurements  from  multiple  source-receiver  pairs.  A  thorough  understanding  of  the  forward 
problem,  including  the  effects  of  horizontal  refraction  and  mode  coupling,  is  necessary  to  successfully 
estimate  environmental  parameters  in  regions  with  3D  inhomogeneities.  This  project  is  funded  by 
ARL:UT’s  IR&D  program. 

Acoustic  propagation  modeling  for  diver  detection  sonar  systems 

The  purpose  of  this  work  is  to  characterize  waveforms  at  virtual  receiver  distances  on  the  order  of  a 
1000  meters  away  from  active  diver  detection  sonar  systems  installed  at  fixed  locations  within 
operational  sites  of  interest.  The  3D  coupled-mode  model  is  applied  for  this  purpose.  This  project  is 
funded  by  Space  and  Naval  Warfare  Systems  Center  Pacific. 
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